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Chapter  1 
Introduction 


The  theoretical  study  presented  in  this  report  has  been  triggered  by  some  experi¬ 
mental  findings  [1]  indicating  a  negative  temperature  ionization  coefficient,  that  is 
a  breakdown  voltage  decreasing  with  temperature.  Traditionally,  semiconductors 
disply  a  positive  temperature  coefficient,  that  is  the  breakdown  bias  increases  with 
temperature.  This  is  a  signature  of  the  intrinsic  stability  of  a  device,  as  under 
high  bias  operation,  the  internal  heating  of  a  device  immediately  translate  into  an 
enhancement  of  the  breakdown  voltage.  On  the  contrary,  a  device  exhibiting  a  neg¬ 
ative  coefficient  is  intrinsically  unstable,  as  an  increase  in  the  internal  temperature 
automatically  favors  the  device  breakdown.  For  instance,  devices  with  a  negative 
temperature  coeffient  would  distructively  react  to  a  sudden  current  glitch.  Further¬ 
more,  the  presence  of  hot  spots,  that  is  special  points  of  a  device  where  (for  instance 
because  of  the  presence  of  dislocations)  a  lot  of  current  is  channeled,  would  also 
have  catastrophic  consequences.  Figure  1.1  shows  a  typical  4H-SiC  PN  structure 
fabricated  by  NASA  using  available  commercial  substrate. 


Figure  1.1:  Cross  section  of  the  jH-SiC  pn  junction  diode 


The  corresponding  IV  characteristic  is  presented  in  Fig.1.2.  In  the  experimental 
investigation  all  diodes  showing  strong  leackage  currents  or  non-vertical  breakdown 
characteristics  have  been  escluded. 
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Figure  1.2:  Curve-tracer  used  for  the  d.c.  IV  characteristics 


The  experimental  technique  was  based  on  the  measurement  of  the  device  response 
to  a  sudden  voltage  pulse  with  a  rise  time  of  about  200  ns.  From  the  temporal 
evolution  of  the  diode  current  and  voltage,  the  temperature  coefficient  has  been 
extracted.  An  sample  of  the  experimental  finding  is  presented  in  fig.1.3.  For  voltage 
valure  below  the  breakdown  voltage,  a  fast  transient  response  is  observed,  which  is 
associated  to  displacement  currents  in  the  device. 

At  higher  voltages,  a  persistent  current  and  voltage  is  found  on  the  device  as 
long  as  the  voltagestep  is  present.  For  4H-SiC  diode  (shown  in  the  figure),  a  positive 
temperature  coefficient  is  found.  On  the  contrary,  all  6H-SiC  sample  that  have  been 
investigated  have  shown  an  opposite  trend. 

The  objective  of  the  present  work  has  been  to  determine  wheter  intrinsic  physical 
phenomena  can  be  responsible  for  the  anomalous  behavior  of  6H-SiC  diodes.  Our 
study  is  based  on  a  selfconsistent  drift-diffusion  approach.  The  physical  basis  of  the 
method  will  be  described  in  Ch.  2,  while  the  numerical  implementation  is  presented 
in  Ch.  3.  The  method  is  then  applied  to  a  SiC  diode  and  the  results  discussed  in 
Ch.  4. 


umo  (n«) 

Figure  1.3:  Experimental  results  on  the  diode  dynamical  response 
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Chapter  2 

Drift-diffusion  approach 


In  this  section,  we  will  describe  the  Drift-diffusion  (DD)  algorithm,  which  constitutes 
the  basis  of  our  simulation. 


2.1  Drift-Diffusion  method 

The  starting  point  for  any  semiclassical  transport  theory  (such  as  the  DD  method) 
is  Boltzmann  transport  equation  which  describes  the  evolution  of  the  carrier  dis¬ 
tribution  function  in  the  presence  of  external  perturbations  and  non  homogeneous 
conditions.  The  importance  of  the  distribution  function  is  that  it  contains  all  mi¬ 
croscopic  information  on  the  carrier  dynamics.  In  particular,  its  first  three  moments 
are  related  to  the  observable  of  the  systems,  namely 

•  carrier  concentration 


”(r>  *)  =  /  /(r>  k>  *)dk  C2-1) 

•  mean  current  density 

J (r>  *)  =  3  /  v(k)/(r,  k,  t)dk  (2.2) 


•  mean  energy 


<E(r’  t))  =  W?  /  E(k)1' (r’ k' t)dk)  (2-3) 

Directly  from  Boltzmann  equation,  it  is  possible  to  obtain  the  differential  equa¬ 
tions  for  such  observables  [2],  namely  the  current  continuity  equation 


dn 

dt 


1  _ 
Vr«7n 

Q 


(2.4) 
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and  the  balance  equation  for  momentum  and  energy, 


djn  =  -2(-g)VrW  g2nF  /  J_\ 
dt  m*  m*  \rm/  n 

dW  /  1  \ 

-w=-ww+w-{-yw-w°)+x 


(2.5) 

(2.6) 


where  W  is  the  kinetic  energy  tensor  (W0  is  its  equilibrium  value),  rm  and  rE  are 
respectively  the  average  time  for  momentum-dissipating  and  energy-dissippating  col¬ 
lisions,  F  is  the  electric  field. 

The  equations  2.4,  2.5  and  2.6  are  the  basis  of  the  so-called  ’’hydrodynamical 
method”  for  device  simulation  (once  coupled  to  Poisson  equation).  A  simpler  model 
can  be  obtained  with  some  assumptions,  as  shown  below.  The  starting  point  is  Eq. 
2.5  expressed  as 


j  1  ~ 2(— g)VrW  fnF 

(1  /An)  dt  m*(l/rm)  m*(l/rm) 


If  we  define  the  mobility  as 

_  q 

^ ”  m*(l/rm ) 

and  assume  that  the  current  does  not  vary  considerably  in  a  time  l/(l/rm),  then 
we  have 


Jn  =  nqnn  F  +  2/znVrW  (2.8) 

The  tensor  W  contains  two  contributions  related  respectively  to  thermal  and  drift 
motion.  In  order  to  simplify  Eq.  2.8  the  drift  component  of  W  are  ignored  and  the 
temperature  tensor  is  assumed  to  be  diagonal,  so  that 


_  tiKbTc  ^ 

¥Vi]  2  —  3 


(2.9) 


where 

W  =  -uKbTc 


(2.10) 


is  the  kinetic  energy  density,  and  Tc  the  carrier  temperature.  By  inserting  Eq.  2.9in 
Eq.  2.8,  we  have 

2 

Jn  =  nqHnF  +  ~/j,nVrW  (2.11) 

By  expanding  the  gradient,  and  expressing  the  result  in  terms  of  Tc,  the  current 
equation  becomes 

Jn  =  nqnn  F  +  qDn'Vn  +  qSnV  Tc  (2.12) 


which  is  indeed  called  ’’Drift-diffusion”  equation.  In  2.12,  Two  new  parameters  have 
been  introduced,  The  diffusion  coefficient: 


Dn 


KbTc 

t^n 

q 


(2.13) 
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and  the  Soret  coefficient: 


(2,14) 


„  _  Kb 
Sn  —  n  fin 
Q 


In  2.12  Tq  must  been  known.  At  low  electric  fields,  the  carrier  temperature  can 
be  taken  equal  to  the  lattice  one;  at  higher  fields,  the  energy  balance  equation  has 
to  be  solved. 


2.2  Generation  terms 

In  Eq.  2.4,  an  extra  term  has  to  be  added  on  the  right  hand  side,  which  account 
for  generation  processes  that  can  alter  the  carrier  concentration.  In  a  device  under 
strong  bias,  two  main  processes  have  to  considered,  impact  ionization  and  field- 
assisted  ionization  via  tunneling. 

2.2.1  Impact  ionization 

In  the  presence  of  high  electric  fields,  carriers  can  gain  sufficient  energy  to  be  able  to 
promote  an  electron  from  the  valence  to  the  conduction  band.  Such  process,  called 
“impact  ionization”,  produce  an  additional  electron- hole  pair  in  the  semiconductor. 
The  multiplication  process,  when  repeated,  can  lead  to  the  breakdown  of  the  device. 

The  rate  of  electron-hole  pair  production  due  to  impact  ionization  processes  G 
is  given  by: 


G  =  annvn  +  appvp 

where  an  is  the  electron  ionization  coefficient  defined  as  the  number  of  electron- 
hole  pairs  generated  per  unit  length  by  an  electron,  vn  is  the  electron  drift  velocity 
and  n  the  electron  concentration.  The  subscript  p  refers  to  holes. 

The  electric  field  dependence  of  the  ionization  coefficients  an>p  can  be  expressed 
as: 


an,p  =  an,pexP 


jpcrit  Q 
(__jhp_Yn’p 

'  jE? 


(2.15) 


where  E  is  the  electric  field,  ECTlt  is  a  critical  electric  field  for  the  ionization  process, 
and  a™p  a  parameter. 


2.2.2  Ionization  due  to  tunneling 

Generally,  in  conventional  semiconductors  at  room  temperature,  dopants  are  fully 
ionized.  In  special  cases,  for  instance  at  very  low  temperature  or  for  dopants  creating 
states  deep  into  the  gap,  ionization  can  be  only  partial.  Then  the  presence  of  an 
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electric  field  can  influence  the  degree  of  ionization  of  dopant  atoms,  via  a  process 
called  field-assisted  tunneling. 

The  electron  generation  rate  due  to  field  assisted  tunneling  can  be  espressed  as: 

^  dnd  nd 

Ort  —  -77—  = - . 

ot  T 

where  nd  is  the  electron  concentration  on  the  donor  level  and  r  is  the  tunneling 
time. 

The  temporal  evolution  of  nd  is  given  simply  by: 

nd(t)  =  ndoe~  &  ^dy 


where  nd0  is  the  initial  concentration  of  electrons  on  the  donor  atoms, 
form: 


ndi  =  ndoie 


In  discrete 


The  electric  field  dependence  of  the  tunneling  time  is  given  by  the  equation[4]: 


where  a  and  ui  are  given  by: 


-1  a  I-2-! 
r  r) 


a 

U) 

Ed  being  the  dopant  energy  level. 


3  eh 
12  Ed 

h  ' 


(2.16) 


2.3  Mobility 


The  field  dependence  of  the  carrier  mobility  /i(F)  is  given  by  the  relation  [5]: 


/Vp(F)  — 


(2.17) 


where  a  is  a  constant  (typically  equal  to  1  or  2).  The  saturation  velocity  vsat  depends 
on  temperature,  according  to  the  phenomenological  law 


Vsat  (T)  = 


Vsat  (0) 

1  +  0.8 exp  ( 600^)  • 


In  SiC  we  have  usot(0)  =  4.64  x  107  cm/s. 


(2.18) 
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>- 


The  temperature  dependence  of  the  carrier  mobility  can  be  expressed  as 


=  umm  + 


..delta 

f  T  \  an'p 

i+wr' 

UooaJ 

(2.19) 


where  the  parameters  p%£a,  A£p,  an,p  e  7„>p  depend  on  the  material. 


2.4  Thermal  effects 

Thermal  effects  due  to  internal  heating  can  be  accounted  for  via  the  heat  transport 
equation:  [6] 

Pc^  -H  =  V-  (k(Tc)VTc)  (2.20) 

where  p  (kg/m3)  and  c  (Joule /kgK)  are  the  density  and  specific  heat,  respectively, 
k(Tc)  the  thermal  conductivity  and  H  the  local  heat,  given  by  the  product  Jn  x  F. 
The  temperature  dependence  of  p  and  c  can  be  neglected  [6].  Furthermore,  if  we  are 
not  interested  in  the  thermal  transient  behavior,  we  can  set  dTc/dt  to  zero,  which 
greatly  facilitate  the  solution  of  Eq.  2.20  in  one  dimension. 


2.5  Displacement  current 

For  the  study  of  electrical  transient,  the  contribution  of  the  displacement  current 
cannot  be  neglected.  For  a  ID  structure  of  length  W,  the  time  dependent  total 
current  density  is 


4  =  qpnnE  +  qDnVn  +  qSnVTc  +  <^L 

W  at 

The  current  can  then  be  written  as 

m  =  c^+ut) 

where  Co  =  is  called  cold  capacitance ,  A  being  the  device  cross  section. 

A  one  dimensional  device  can  therefore  be  represented  as  a  current  generator  in 
parallel  to  the  cold  capacitance.  At  every  time  step  in  the  simulation,  the  voltage 
V  must  be  updated,  as  it  depends  on  the  charge  stored  in  the  capacitance. 
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Chapter  3 

Numerical  model 


In  this  chapter  we  will  outline  the  physical  approach  used  in  the  study  and  its 
numerical  implementation 

3.1  Fundamental  equations 

As  pointed  out  in  the  previous  section,  the  fundamental  equations  of  our  Drift- 
Diffusion  approach  are  [6] 

•  Poisson  equation 

V-(eV$)  =  -q(p-n  +  N%  -7VJ)  (3.1) 

•  continuity  equation  for  electrons  and  holes 

_  T  <9n  _ 

v'-7"-«¥=?b 

V  '  +  q~  ~  —qR  (3.2) 

•  current  equation  for  electrons  and  holes 

Jn  =  nqpnF  +  qDnVn  +  qSnVTc 
Jp  =  pqiJ,pF  -  qDpVp  -  qSpVTc  (3.3) 

•  heat  transport  equation 

dTc 

pc- -  J  x  F  =  V  •  0 k(Tc)VTc )  (3.4) 


The  models  for  R,  the  generation-recombination  rate,  for  pn,pp  (electron  and  hole 
mobility)  and  for  K(T)  will  be  discussed  below. 


9 


3.2  Discretization 

Here  we  will  discuss  the  various  discretization  techniques  followed. 


3.2.1  Continuity  equation 

By  substituting  Eqs.  3.3  in  3.2,  the  following  relations  are  found 


dn(x,t )  _  .  ,dn(x,t) 

—  Vn  (£,  t)  - 


dt 


dx 


dx 


dx 2 


dn(x,  t)  dDn(x,t)  _  c  d2Tc(x,t)  t  dTc(x)  dSn(x) 

*  & n  rx  o 


dx 


dx 


dx2 


dx  dx 


(3.5) 


MM  =  ^^MM  t)  +£)/p(M 


dt 

+ 


<9x 


dx 


dx2 


dpfa  t)  dDpjx,  t)  +  s  &Tc{x,t)  ___  d7c(x)  dffpfx) 


dx 


dx 


dx2 


+ 


dx  dx 


(3.6) 


where  vnjP(x)  are  the  carrier  velocities  defined  as: 

Vn  =  “//n  F 


Vp  =  /ipF 


(3.7) 

(3.8) 


For  the  mobility  we  have  adopted  either  of  the  two  models  decribed  earlier.  In 
particular,  we  have  used  the  model  2.17  for  fixed  temperature  simulations,  and  that 
of  2.19  when  the  time  and/or  space  temperature  variation  has  been  accounted  for. 
The  discetization  of  3.5  e  3.6  has  been  obtained  with  the  following  approximations 
for  the  spacial  derivatives: 

fi  =  12A2  ~~  +  8^i+1  ~~ 

fi  ~  i2Ax2  +  ~  30/i  +  16/i+i  -  /*+: 2) 


while  for  the  time  derivative  a  simpler  discretization  has  been  used: 


— -fo+j  +  ni-i) 
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3.2.2  Poisson  equation 

During  the  simulation,  it  is  necessary  to  upgrade  the  electric  field  value  at  fixed 
times,  whose  separation  is  determined  by  physical  consideration  (e.g.  the  plasma 
frequency  of  the  carriers). 

In  one  dimension,  Poisson  equation  takes  the  form 

(3.9) 

where  e  is  the  semiconductor  dielectric  constant  and  $  the  electrostatic  potential, 
while  the  net  charge  p  is  given  by 


p(x)  =  -q[n(x)  - p(x)  -  N£  +  Na]  (3.10) 

where  q  is  the  electron  charge,  iVj  and  are  respectively  the  ionized  acceptor  and 
donor  densities,  and  n  e  p  are  the  free  carrier  (electron  and  hole)  concentrations. 
The  ionized  donor  (N£)  and  acceptor  (iVj)  concentrations  are  given  by 


nd  = 

Na  = 


Nd 

l  +  gdexp(2j=£2) 
Na 

1  +  gvexp(  K~bTf ) 


(3.11) 


where  Ef  is  the  Fermi  level,  Ed  and  Ea  the  donor  and  acceptor  energy,  and  =  2, 
gv  =  4  their  degeneracy  factors. 


Among  the  several  ways  to  solve  Poisson  equation  we  have  adopted  the  so-called 
’’Gauss  elimination”  technique.  We  subdivide  our  ID  system  into  N  cells  (Poisson 
lattice),  and  discretize  the  equation: 

$i-i  -2$i  +  $i+1  =  -Ax2^  (3.12) 


where  Ax  is  the  fixed  dimension  of  each  cell.  Boundary  conditions  are  imposed  on 
two  extra  cells,  as  shown  in  fig.  3.1 

$o  = 

$N+ 1  —  Vfin 

Vin  and  Vfin  being  the  bias  levels.  Let  us  consider  the  generalization  of  Eq.  (3.12) 

+  o$i  +  a+$j+j  =  p^  (3.13) 


and  set 

$t+l  =  (*»$»  +  Pi 


(3.14) 
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where  a*  and  /%  are  constants.  By  substituting  the  previous  expression  into  (3.13) 
we  obtain: 


—  0>  1  ”b  liifl+Pi  Pi)  (3.15) 


with 


7 i  = 


1 

a  + 


(3.16) 


By  comparing  Eq.  (3.15)  with  Eq.  (3.13)  the  expression  for  the  constant  coefficients 
is  found: 


&i— 1  — 

Pi-i  =  7 i(a+Pi-p'i)  (3.17) 

By  imposing  the  first  boundary  condition 

=  V] fin  —  +  fix  (3.18) 

we  have 

o;jv  =  0 

Pn  —  V fin-  (3.19) 

The  method  is  now  clear:  starting  from  the  condition  (3.19)  and  using  Eqs.  (3.17) 
the  on,  /3i  coefficients  are  determined.  Using  then  the  second  boundary  condition 
($o  =  Vin)  Eq.  (3.14)  gives  the  potential  values  over  the  whole  lattice. 

In  the  special  case  of  interest,  that  is  the  solution  of  Poisson  equation: 

a_  =  a+  =  1 


a  =  — 2 
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On  scalar  computers,  the  Gauss  elimination  method  is  much  faster  than  matrix 
inversion  approaches,  as  it  requires  N  operation  instead  of  N2. 

Once  the  electrostatic  potential  profile  is  calculated,  the  electric  field  is  found  from 

(3.20) 

which  in  the  discrete  version  becomes: 


Ft  = 


$t-2  ~  8$i-l  +  8$M-i  —  2 

12Ax 


(3.21) 


In  order  for  the  solution  to  have  physical  meaning,  the  parameters  adopted  in 
the  numerical  solution  of  any  differential  equation  have  to  obey  some  constrains. 
For  our  case  the  mesh  size  has  to  be  lower  then  the  Debye  length,  given  by: 


where  N  is  the  doping. 


Ld  = 


{ 


ckBT 

q*N 


by 


The  temporal  time  step  has  to  be  smaller  of  the  dielectric  relaxation  time,  given 


tdr  — 


qNfi' 


3.2.3  Heat  transport 

The  equation  describing  heat  transpot  (Eq.  3.4): 

pc^  -  J  x  F  =  V  •  (J k(Tc)VTc )  (3.22) 


can  be  discretized  as  follows  (assuming  a  constant  thermal  conductivity  k): 


prn+l  __  pm  ^ 


pc 


-Til*  +  16(^-im  +  T&y)  -  307?”  -  7^2  J?*F? 

12Ax(i)  +  k 


At 

(i  =  1,..,N), 

from  which  the  temperature  in  the  i-th  point  at  the  (m  +  l)-th  instant  is  found: 
k*  At 


pTTl+1  =  rpm  _|_ 


pc 


-T£2  +  16(772!  +  TftJ  -  307)m  -  T£2  J T  *  F? 

12A®(»)  +  k 


(i  =  l,..,N), 
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where  Jf1  and  F™  are  the  total  current  density  and  the  electric  field  in  the  i-th 
point.  The  boundary  conditions  are 

T(-l)  =  T(0)  =  T(N  +  1)  =  T(N  +  2)  =  300  K 


3.3  Generation  processes 

Two  generation  processes  are  important  in  our  study:  impact  ionization  and  tun¬ 
neling.  As  for  the  former  process,  in  the  absence  of  experimental  values  for  the 
4H  politype,  we  will  take  the  same  values  as  for  the  6H  one.  This  assumption  is 
justified  by  the  fact  that  the  breakdown  voltages  for  PIN  diodes  obtained  from  the 
two  materials  are  quite  similar. 


For  the  a-SiC,  at  room  temperature  impact  ionization  is  caused  preferably  by 
holes.  Since  the  valence  band  structure  is  similar[7],  the  same  should  occur  in  /?-SiC. 
At  higher  temperatures,  electrons  also  start  contributing  to  the  ionization,  and  a 
different  behavior  of  the  two  politype  has  to  be  expected  due  to  the  differences  in 
their  conduction  bands  [7]. 

The  total  number  of  electron-hole  pairs  generated  per  unit  time  is 

Gn  =  ~{&nJn  "t"  OipJp ) 

where  an  and  ap  are  tensors,  with 

an||  =  an,  CKP||  =  ap 


and 


otn±  —  <an/3.5, 


=  0£p 


otn>P{E)  =  an,pexp  -  j  . 

On  the  basis  of  a  fit  with  available  experimental  data,  an  and  ap  have  the 
following  dependence  on  temperature: 

an(T)  =  7.628  x  106(T/300)2  +  1.506  x  106(T/300)  -  4.184  x  lO^cm"1] 

aP(T )  =  1.469  x  106(T/300)2  +  1.083  x  107(T/300)  -  3.100  x  lO^cm"1] 

bn(T)  =  2.58  x  106  [V/cm]  (indipendent  of  T) 
bp(T )  =  1.90  x  107  [V/cm]  (indipendent  of  T) [7] . 


As  for  the  tunneling  contribution  ( Gn  and  Gn  for  electrons  and  holes,  respec¬ 
tively)  we  have  used  the  following  model: 

f~i  _ 

Un,P  ~ 

'n,p 
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3.4  Basic  structure  of  the  simulated  N+NP+  diode 

The  structure  of  the  simulated  diode  is  shown  in  Fig.  3.2.  The  diode  is  made  up  by 
a  0.25  pm  wide  N-region  (doped  at  518cm-3),  followed  by  a  second  N-type  region  of 
0.6  fj,m  doped  at  1018  cm-3. 

N  + 

N 

P  + 

Figure  3.2:  Simulated  diode  structure 

The  third  region  is  p-type,  with  a  density  of  520,cm-3  and  a  width  of  0.3  /xm. 
The  contacts  on  the  heavily  doped  n  and  p  regiona  are  made  of  Aluminum. 
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The  dimension  and  doping  levels  of  the  device  have  strong  influence  on  the 
simulation.  In  fact,  in  order  to  have  a  sufficient  number  of  mesh  points  a  spacial 
mesh  Ax  =  10_9m  has  been  taken,  which  has  in  turn  implied  a  large  number  of 
discretization  points  and  a  heavy  computational  load. 

Furthermore,  in  order  to  guarantee  the  stability  of  the  finite  difference  algorithm, 
the  Von  Neumann  criterion  requires  that 


vAt 

Ax 


<  1. 


Since  the  saturation  velocity  in  SiC  is  2  x  107cm/s,  a  value  At  <  5  x  10-14s  is 
needed.  Actually,  an  even  smaller  time  step  is  requires  by  the  fact  that  tunneling 
times  as  short  as  5  x  10-14  are  calculated  in  the  high  field  region  of  the  device,  which 
forced  us  to  use  time  steps  of  the  order  of  2  x  10-16,  thus  increasing  enormously  the 
complexity  of  the  numerical  algorithm.  In  order  to  speed  up  the  simulation  under 

transient  conditions,  the  depletion  approximation  for  each  junction  is  used  as  initial 
condition  to  each  simulation. 


3.5  Simulation  parameters 

The  following  table  (tab.  3.1)  illustrates  the  parameters  used  in  the  simulation  of 
6H  e  4H-SiC  diodes. 
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Table  3.1:  Simulation  parameters 


6  H 

AH 

Dielet. const. e 

9.66 

9.66 

Energy  gap 

!S£ng| 

Electron  mass 

mn0  =  0.41 

3 

3 

o 

II 

O 

t— i 

Hole  mass 

mp  o  =  1 

rripo  =  1 

Acceptor  level 

Ea  =  0.2eV 

Ea  =  0.191eV 

Donor  level 

Ed  =  O.leV 

Ed  =  0.065eV 

Ioniz.  const. 

CM  Th 

II  II 

9c  =  2 

=  4 

(cm2  /Vs) 

nn 

(cm~3) 

RjBQQI 

HBUbB 

9 

6  H 

4# 

QLn  —  3 

Oip  —  3 

O^n  =  3 

Oip  =  3 

7n  =  0.45 

7?  =  0.5 

7n  =  0.45 

7P  =  0.5 

Pn  =  2 

/?n  =  2 

Sat.  vel.  (cm/s) 

■HOyyftljH 

K2SIIBB 

Impact  ionization 

(cm  *) 

an  =  1.66  x  106 
ap  =  5.18  x  106 

a„  =  1.66  x  106 
op  =  5.18  x  106 

(V/cm) 

wM 

1 

Thermal  cond. 
W/cmK 

AT  =  4.9 

a:  =  4.9 
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Chapter  4 

Diode  simulation 


As  a  first  step  for  the  simulation,  we  have  verified  that  the  models  provides  realistic 
profiles  under  ’’normal”  conditions.  For  instance,  we  have  biased  the  diode  under 
reverse  voltage  and  checked  the  with  of  the  depletion  region  with  respect  to  analitical 
results.  Once  the  physical  model  and  the  numerical  algorithm  have  proven  to  be 
correct,  we  have  started  a  series  of  simulation  of  the  time  dependent  response  of  SiC 
diodes  to  a  voltage  step.  Due  to  the  heavy  computational  load,  we  have  not  been 
able  to  apply  rising  times  and  pulse  of  ns  duration,  as  in  the  experiment  of  Fazi  and 
Neudeck.  Rather,  we  have  simulated  rise  times  of  the  order  of  a  ps,  and  we  have 
looked  at  the  diode  dynamics  toward  the  steady  state  condition. 

4.0.1  Diode  response  in  the  absence  of  field-assisted  ioniza¬ 
tion  and  impact  ionization 

In  figs.  4.1(a)  and  (b)  we  show  the  time  evolution  of  the  applied  voltage  {Vb,in  —  Vapp) 
and  of  the  electron  current  density  Jn.  The  application  of  a  reverse  voltage  step 
of  100V  over  1  ps  (from  t=4ps  up  to  t=5ps)generates  a  current  pulse  due  to  the 
carrier  movement  following  the  sudden  change  of  the  electric  field.  At  the  end  of 
such  transient,  the  current  goes  back  to  zero,  as  expected  under  steady  reverse  bias 
in  the  presence  only  of  thermal  effects. 

In  fig.4.2  we  see  that  the  electric  field  evolves  towards  a  steady  state  value  much 
higher  than  that  reached  at  the  end  of  the  voltage  step.  Thus,  the  transient  does 
not  lead  to  an  overshoot  in  the  field  peak  value,  as  it  might  have  been  expected. 
The  dynamics  of  the  carrier  concentratio  (Fig.  4.3)  is  consistent  with  the  electric 
field  one,  showing  the  large  depletion  region  reached  in  steady  state. 

4.0.2  Diode  response  in  the  presence  of  tunneling  or  impact 
ionization 

When  tunneling  is  considered,  the  current  response  has  the  same  shape  as  seen  before 
(Fig.4.4),  but  the  peak  value  is  much  higher.  This  is  due  to  the  fact  that  the  high 
field  values  reached  while  approaching  stationary  conditions  lead  to  sizeable  carrier 
tunneling  out  of  donor  and  acceptor  states.  This  is  confirmed  by  the  distributions 
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Electron  &  Hole  density  (cm'3)  Field  (V/cm) 


Figure  4.2:  Electric  field  evolution 


Figure  4.3:  Carrier  evolution  during  transient 
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of  carriers,  unionized  donors  and  acceptors  (fig.4.7).  At  the  end  of  the  voltage  step, 
all  donors  in  the  N  region  are  completely  ionized,  while  in  the  N+  and  P+  the 
complete  ionization  condition  is  not  reached.  The  net  result  is  that  the  electron  and 
hole  concentration  are  respectively  a  factor  two  and  three  than  in  the  absence  of 
tunneling. 


6e-i-06 


18.0 


Figure  4.4:  Current  response  in  the  presence  of  tunneling 


As  a  consequence  of  the  increased  amount  of  free  and  fixed  charges,  the  elec¬ 
tric  field  values  are  also  higher  (the  maximum  field  grows  from  3  x  106  V/cm  to 
5.5  x  106  V/cm).  From  the  time  evolution  of  the  various  quantities,  it  can  be  no¬ 
ticed  that  the  tunneling  process  is  very  fast,  as  espected  from  the  large  value  of  the 
corresponding  rates. 

i,From  Fig.  4.7  we  notice  the  complete  donor  ionization  in  the  n  layer,  and  the 
only  partial  ionization  of  the  dopants  in  the  highly  doped  regions. 

The  electron  density  profile  at  the  N+N  displays  at  very  short  time  an  interesting 
feature.  As  the  electrons  are  injected  in  the  N-zone,  they  experience  a  strong  ac¬ 
celeration  towards  the  N+  region  caused  by  the  strong  electric  field  of  the  depleted 
zone.  As  they  reach  the  interface,  the  electric  field  abruptly  drops  and  the  electrons 
slow  down,  leading  to  the  accumulation  layer  visible  in  the  picture.  In  order  to  check 
if  such  short-time  effect  still  exists  when  the  rise  time  of  the  voltage  step  is  longer, 
we  have  simulated  a  single  voltage  step  of  100F  with  a  rise  time  of  100  ps  starting 
from  t=0. 


21 


Figure  4.5:  Electric  field  evolution  in  the  presence  of  tunneling 


Figure  4.6:  Carrier  concentration  in  the  presence  of  tunneling 


Accettori  non  ionizzati  (cm 3) 


Figure  4.7:  Unionized  donor  (a)  and  acceptor  (b)  concentration  in  the  presence  of 
tunneling 


Figs. 4. 8, 4. 9,4. 10  show  the  temporal  evolution  of  the  current,  carrier  concentra¬ 
tion,  and  electric  field,  respectively,  indicating  the  the  electron  acculation  found 
previously  was  due  to  the  ultrashort  rise  time.  In  steady  state,  very  similar  results 
are  found. 


Figure  4.8:  Current  evolution  for  a  rise  time  trise  =  lOOps  (tunneling  only) 


Figure  4.9:  Carrier  concentration  evolution  for  trise  =  lOOps  (tunneling  only) 
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Figure  4.10:  Electric  field  evolution  fortrise  =  lOOps  (tunneling  only) 
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When  only  impact  ionization  is  considered,  no  appreciable  different  is  found  with 
respect  to  the  case  of  no  ionization.  That  is,  the  electric  field  is  always  below  the 
threshold  value  for  such  process.  As  a  result,  the  current  goes  to  zero  after  the  initial 
transient  (fig.4.11). 


3.0  8.0  13.0  18.0 

time  (ps) 


Figure  4.11:  Corrente  in  the  presence  of  impact  ionization  (Vun  ~  Vapp  —  100 W ) 
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12:  Electric  field  in  the  presence  of  impact  ionization  (Vb.in  —  Vapp  = 


Figure  4.13:  Carrier  concentration  in  the  presence  of  impact  ionization  (Vb. 


Figure  4.14  shows  that  a  non  vanish  (starting  at  T=4ps)  steady  current  is  ob¬ 
tained  when  a  bias  of  170  V  is  applied  with  a  rise  time  of  2  ps.  Thus,  the  contribution 
of  impact  ionization  to  the  persistent  current  is  actually  a  function  of  the  value  and 
slope  of  the  initial  potential  step. 


time  (ps) 


Figure  4.14:  Current  density  in  the  presence  of  impact  ionization  (V^m  —  Vapp  = 
170  V ) 


Dynamically,  we  can  see  that  while  the  current  tends  to  return  to  zero,  ionization 
processes  start  that  push  the  current  value  up  to  its  stationary  value.  The  profiles  of 
carrier  concentration  (figs. 4. 16  and  4.17  indicate  that  significant  ionization  (and  the 
corresponding  widening  of  the  depletion  region)  initiates  around  1  ps  after  the  step. 
The  electric  field  displays  an  enhancement  in  the  N  region  and  a  decrease  in  both 
N+  e  P+  zones  during  the  initial  current  pulse,  followed  by  a  reversed  behavior  when 
impact  ionization  start  contributing.  From  figs.4.18  (a,b)  we  can  see  that  indeed 
the  ionization  rate  achieves  its  maximum  value  after  12ps. 
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Figure  4.15:  Electric  field  in  the  presence  of  impact  ionization  (Vb,in 
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Figure  4.16:  Electron  concentration  in  the  presence  of  impact  ionization  (Vbin 
Vapp  =  170V) 
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Figure  4.17:  Hole  concentration  in  the  presence  of  impact  ionization  fHm  ~  Va 
170  V) 
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Figure  4.18:  Impact  ionization  rate  (Vb. 


4.0.3  Response  in  the  presence  of  tunneling  and  impact  ion¬ 
ization 

When  both  tunneling  and  ionization  are  present,  a  finite  stationary  current  is  found, 
in  agreement  with  the  experimental  findings.  The  current  shows  an  oscillatory  be¬ 
havior  in  response  to  the  step  potential,  but  does  not  return  to  zero  as  it  is  substained 
by  a  continious  generation  of  carriers  which  eventually  leads  to  the  disappearence 
of  the  depletion  region. 


5e+06 


2.0 


12.0 


22.0 
time  (ps) 


32.0 


42.0 


Figure  4.19:  Current  density  in  the  presence  of  tunneling  and  ionization 


At  the  end  of  the  pulse,  the  ionization  rate  is  fairly  low,  but  it  grows  steadily 
with  time  as  the  electric  field  increases.  In  the  proximity  of  the  junction  the  field  is 
actually  so  high  to  give  almost  complete  ionization  of  the  dopants. 

In  figs.4.19,4.20,4.21,4.23  the  temporal  evolution  of  electric  field,  carrier  concen¬ 
tration,  and  depletion  charge  is  reported.  Near  the  N  -  P+  the  hole  concentration 
presents  a  peak  due  to  the  large  contribution  of  impact  ionization  processes  caused 
by  the  junction  field. 
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Depletion  charge  (cm'3)  Depletion  charge  (cm'3) 


Figure  4.23:  Fixed  charge  in  the  presence  of  tunneling  and  ionization 
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4.0.4  Influence  of  deep  levels 

We  have  also  studied  the  influence  of  deep  levels  on  the  electrical  response  of  SiC 
N+NP+  diode.  The  doping  levels  and  the  width  of  each  region  are  respectively 
5  x  1018,  1018,  5  x  102°  cm-3  and  0.2,  0.6  e  0.3 \xm.  A  rise  time  of  lOOps  has  been 
considered. 

Figure  4.24  (dashed  line  shows  the  current  behavior  for  donor  and  acceptor  levels 
respectively  at  Ed  =  100  meV,  Ea  =  600  meV.  Only  impact  ionization  processes 
are  considered.  The  initial  current  spike  is  due  to  a  displacement  contribution.  The 
current  then  raises  further  as  the  internal  electric  field  grows  towards  the  cathod, 
leading  to  a  rather  weak  generation  of  holes.  Most  acceptors  remain  unionized. 
The  slow  current  rise  reflects  well  the  experimental  findings,  indicating  that  deep 
acceptor  levels  do  indeed  contribute  to  the  conduction  mechanisms. 


Figure  4.24:  Electrical  current  response  in  the  presence  of  deep  levels 


When  tunneling  is  also  considered,  the  current  does  not  tend  to  saturate  but 
inreases  further,  up  to  values  a  factor  three  higher  than  in  the  absence  of  tunneling 
processes.  Such  enhancement  is  directly  related  to  the  increased  number  of  free 
carriers  and  ionized  acceptors  (see  fig.4.24)  in  the  device,  favoured  by  the  high 
electric  field  at  the  edge  of  the  N+  and  P+  zones. 

4.0.5  Internal  heating  effects 

We  have  perviously  described  the  simplified  model  that  allows  us  to  account  for 
the  temperature  variation  connected  to  the  current  flux  inside  the  diode,  and  to 
calculate  its  profile.  The  algorithm  that  we  have  implementented  gives  us  directly 
the  temperature  map  inside  the  device  at  every  time  step,  thus  providing  information 
on  the  influence  of  internal  heating  on  the  dynamical  response  to  the  voltage  pulses. 
A  temperature-dependent  mobility  model  is  used  to  update  constantly  the  carrier 
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Figure  4.25: 
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Figure  4.26:  Hole  concentration 
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Figure  4.27:  Unionized  donor  concentration 
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Figure  4.28:  Unionized  acceptor  concentration 


mobility  to  the  changing  temperature.  In  the  simulation  a  voltage  step  of  100  V  has 
been  considered  with  a  rise  time  of  lps (star ting  at  4ps). 


Figure  4.29:  Electrical  current  in  the  presence  of  self  heating 


Figures  4.29  and  4.30  show  respectively  the  evolution  of  the  electrical  current 
and  of  the  internal  temperature.  Clearly,  the  temperature  reacts  much  slower  than 
the  current  to  the  voltage  pulse,  a  consequence  of  the  longer  characteristic  times  for 
heat  conduction. 

The  effect  of  the  increasing  temperature  can  be  seen  from  figs.4.32  and  figg.4.33 
e  4.34,  showing  respectively  the  carrier  concentration  and  drift  velocity,  for  both 
electrons  and  holes.  The  temperature  increase  causes  a  reduction  in  time  of  the 
carrier  mobility  and  drift  velocity.  The  carrier  concentration  in  the  doped  contact 
regions  increases  due  to  the  higher  thermal  ionization  of  the  dopants. 
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Figure  4.30:  Time  evolution  of  the  internal  temperature  during  the  current  tr 


Figure  4.31:  Electric  field  evolution 
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Figure  4.32:  (a)  Transient  evolution  of  the  carrier  concentration;  (b)  Approach  to 
the  stationary  condition 
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Figure  4.33:  Electron  velocity  evolution 


4.0.6  Simulation  of  a  l^H-SiC  diode 

We  have  performed  a  simulation  of  a  J^H-SiC  diode,  in  order  to  check  if  there  are 
significant  differences  between  the  two  polytype.  The  electrical  response  of  the  two 
system  has  been  found  to  be  very  similar  (fig.  4.35).  With  respect  to  the  6H, 
the  4H-SiC  displays  a  larger  tunneling  contribution  due  to  the  fact  that  the  donor 
and  acceptor  levels  are  shallower  (Ed  =  0.065eV  and  Ea  —  0.191eV).  The  carrier 
mobility  of  the  two  material  is  also  different: 

/r„o  =  360  err? fV sec  and  [in0  =  800  cm2 /V sec 

lipQ  =  60  cm2 /V sec  e  fj,po  =  100  cm2 /V sec 
.  for  6H-SiC  and  4H-SiC,  respectively. 
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Figure  4.34:  Hole  velocity  evolution 


4.1  DC  I-V  characteristics 


The  stationary  current-voltage  characteristics  for  6H-SiC  extracted  from  the  time 
dependent  simulations  are  presented  in  fig.4.36). 


Figure  4.36:  I-V  Characetristics 


The  two  curves  refer  to  a  simulation  at  constant  temperature  (solid  line)  and 
one  where  internal  heating  is  accounted  for  (dashed  line).  While  the  onset  of  the 
breakdown  is  similar  in  the  two  cases,  when  the  internal  temperature  is  considered 
we  observe  a  slower  current  rise. 

The  actual  value  of  the  breakdown  voltage  is  lower  than  the  experimental  one 
due  to  the  smaller  dimension  of  the  diode  considered  in  the  simulation. 

From  fig.  4.37  we  can  infer  that  at  higher  temperatures  the  breakdown  voltage 
tens  to  incraese.  Thus,  the  conlcusion  of  our  study  is  that  the  negative  tempera¬ 
ture  inization  coeffiecient  resulting  from  the  experimental  analysis  of  6H-SiC  diodes 
cannot  be  attributed  to  internal  physical  mechanisms,  but  probably  to  technological 
problems. 
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Figure  4.37:  IV  characteristics  of  a  6H-SiC  diode  at  different  temperatures 
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